Peierls instability, periodic Bose-Einstein condensates and density waves in 
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We study the quasi-one-dimensional (Q1D) spin-polarized bose-fermi mixture of atomic gases at 
zero temperature. Bosonic excitation spectra are calculated in random phase approximation on the 
ground state with the uniform BEC, and the Peierls instabilities are shown to appear in bosonic 
collective excitation modes with wave- number 2k f by the coupling between the Bogoliubov-phonon 
mode of bosonic atoms and the fermion particle-hole excitations. The ground-state properties are 
calculated in the variational method, and, corresponding to the Peierls instability, the state with a 
periodic BEC and fermionic density waves with the period n/kF are shown to have a lower energy 
than the uniform one. We also briefly discuss the Q1D system confined in a harmonic oscillator 
(HO) potential and derive the Peierls instability condition for it. 

PACS numbers: 03.75.Fi, 05.30.Fk,67.60.-g 
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I. INTRODUCTION 

The cooling and trappingtechniques of atoms, used 
in the realization of BEC Q], and degenerate fermi 
gases |3(, have recently produced quantum bose-fermi 
mixtures of dilute atomic gases 0, ||- Motivated by 
the controllability of various experimental parameters 
such as trap geometries, particle numbers, and interac- 
tion strengths in experiments of the ultra-cold atoms, a 
number of theoretical studies have been done to explain 
various ex peri mental results done for bose-fermi mixtures 
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Further experimental advances have led to the real- 
ization of BEC in quasi-one-dimensional (Q1D) system 
H d, E|> by putting atoms in cylindrical traps long- 
enough that the one-particle energy-level spacing in ra- 
dial direction exceeds the interatomic-interaction energy, 
and the atoms can move effectively in the axial direction. 
Characteristic properties of the Q1D BEC have been 
demonstrated in the formation experiments of matter- 
wave solitons in repulsively/ attractively interacting sys- 
tems [20L l2l| . Also several new phenomena, which are 
typical in the Q1D bose system, have been proposed the- 
oretically: quasi-condensates with fluctuating phases |22j 
and a Tonks- Girardeau gas of impenetrable bosons [23 
and so on. 

In electric conductors, the Q1D systems have been re- 
alized in systems where electric currents can flow easier 
in one s peci fic crystal direction than in the vertical di- 
rections [24j. One of the most fascinating phenomena in 
such a system is the "Peierls instability" : a lattice dis- 
tortion with the period of 2Uf This instability gives 
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rise to a charge-density wave in the ground state, a col- 
lective state of electrons, which is caused by the strong 
correlations among electrons due to the phonon-electron 
interaction. The origin of the effect is the divergence 
in phonon spectrum by the low-energy particle-hole (p- 
h) excitations near fermi surface with the wave-number 
2Uf j which give large contribution in phonon self-energy 
because of a phase-space reduction in Q1D systems. 

In the present paper, we investigate the occurrence 
of the Peierls instability in the Q1D uniform system of 
atomic-gas bose-fermi mixture at T = 0. Different from 
conductors that have the periodic structures by lattices, 
it have no periodic structures originally, and the insta- 
bility and the periodic structure in the ground state is 
brought purely by the coupling between the Bogoliubov- 
phonon mode of bosonic atoms and fermion p-h excita- 
tions. 



This paper is organized as follows. In Section II, we 
introduce the model and derive the equations in random 
phase approximation using the Green's function method 
to calculate the bosonic excitation spectra. In Section 
III, the equations obtained in Section II are applied for 
the system with the uniform ground state, and solved 
numerically. We show the calculated bosonic excitation 
spectra and that the Peierls instability appears in collec- 
tive excitation modes around the wave-number 2fcj?. The 
analytical estimations are also given for the collective ex- 
citation spectra. In Section IV, we show that the state 
with periodic BEC and the density-wave states for both 
boson and fermion (with the period of n/kp) is more 
stable than the uniform state. In Section V, we briefly 
discuss the Peierls instability in the finite Q1D system 
confined in a HO potential 
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II. MODEL, GREEN'S FUNCTIONS AND 
MEAN-FIELD BASIS 

We consider a system of spin-polarized atoms at T = 0: 
N b bosons and Nf fermions with masses m b j each other, 
confined in an axially-symmetric HO potential U b j = 
\m b j {lj 2 {x 2 + y 2 ) + u>lz 2 } where oj r > uj a . When 
atoms are confined tightly enough in the r-direction that 
the radial (r) part of their wave functions can be approx- 
imated by that of the lowest-energy state in the 2D HO 
potential ^m,bjco 2 (x 2 + y 2 ), then we call that the system 
is in the Q1D region. It should be realized when, 

• The r-HO quanta fiu r should be higher than the 
boson-boson/boson- fermion interaction energies: 
ejj^ = 4nh 2 a bb nl D /nib and e b f* = 2Trh 2 abfn 3 c D /m r , 
where n b D j are the boson/fermion 3D-number den- 
sities and abb,bf are the boson-boson/boson- fermion 
s-wave scattering lengths. The m r in e b f is a re- 
duced mass: m r = mbirif /(m,b + to/). 

• The r-HO quanta should also be higher 
than the fermi energy obtained by e|r = 
h 2 {6TT 2 n 3 f D ) 2 / 3 /2m f 

The atomic states for the axial (a) degree of free- 
dom are described by the z-dependent field operators for 
bosons and fermions, <j>{z) and ip(z); the effective Hamil- 
tonian for them is modeled by 



H = Jdz ${z) 
dz ?/v (z 



h 1 d 2 

2to/, dz 2 



+ Ub(z) - fib 



h 2 d 2 

2m f dz 2 



+ Uf(z)-(i f 



dz (fy{z) 



g f^(z)4>(z) 



+9bf4>\z)4>(z) 



Hz), 



(1) 



where the Ubj are the ID chemical potentials for bosons 
and fermions. The effective ID coupling constants gbb 
(boson-boson) and g b f (boson- fermion) in l|14f> are related 
to the s-wave scattering lengths: gbb — 2hui r a bb , 9b f = 
2hu> r abf, which can be obtained from the 3D interaction 
by integrating out the radial parts of the wave functions 
[T3, [23 • It should be noted that, in the limit w a — > 0, we 
obtain the uniform system in the axial dimension. 

In weak interacting systems, the mean field approxi- 
mation should give a good starting point for further cal- 
culations, where bosons are assumed to occupy the lowest 
single-particle state <fo(z) at T = 0, and the expectation 
value of the boson field-operator (order parameter for 
BEC), $(z) = (4>(z)} = ^ip {z), satisfies ID Gross- 
Pitaevskii equation: 



ti 2 d 2 
2mb dz 2 



U b (z) + g bb n b + g b fUf 



$ = /i b $, (2) 



where the (axial) ID boson density n b is defined by 
n b (z) = \<P(z)\ 2 . 

For fermions, the ground state is obtained by a Slater 
determinant of fermion single-particle wave functions 
ipj(z) for the mean-field energies ej below the ID fermi 
energy ep — [if, which are obtained from the Hartree 
equation: 



2to / dz 2 



Uf(z) + g b fn b (z) 



^(z)=e f Mz). (3) 



In this state, the fermion density is obtained by nj(z) — 

E £ /< £F I^-I 2 - 

j 

We are interested in the bosonic excitation modes of 
the system, which can be obtained by the bosonic Green's 
functions 



iG{xi,x 2 ) 

( (T[<p( Xl )&(x2)]) (T[0( X1 )0(X2)]) 
[ (T[^( Xl )^(x 2 )}) (T[tf( Xl )<p{x 2 )]) 



(4) 



where x — [z,t) and <p(x) = <fi( z ,t) — $(2). They satisfy 
the Dyson equation: 



d 



ih—T 3 -K°{ Zl )l 



G(xi,x 2 ) = hS(xi — x 2 )l 
J d 2 x 3 k'E(xi : x 3 ) ■ G(x 3 ,x 2 ), (5) 



where K b (z) = — {Ti 2 /2m b )d 2 / dz 2 + U b {z)~ fi b is a single- 
particle operator and the matrices T3 and 1 are defined 
by 



T3 











u -1)' 


1 = 


i 





(6) 



The boson self-energy in Eq. JSJ is in the 2 x 2-matrix 
form, and separated into the static and the dynamical 
parts: S(a;i, £2) = So(xi, £2) + H(xi, £2). In the mean- 
field approximation, the static part is given by 



, 2n b + Zy-rif 
h-£ =g bb [ 3 



2nb + gin, ) ^ - X2) ' 



In the dynamical part, we take the contributions from the 
fermion particle- hole (p-h) pair-excitations (Figure. ^): 



Kn.{xi,x 2 ) = - 



$(zi)$*(z 2 ) $(zi)$(z 2 ) 



h 9b f y $*( Zl )$*( Z2 ) $*( Zl )$( Z2 ) 
xG- f (x 1 ,x 2 )G^a: 2 ,a ;i ), (8) 



which is just a polarization potential for pair-excitations. 
It should be noted that this approximation is equiva- 
lent to the random phase approximation, and gives the 
strong correlation effects, especially in the bosonic col- 
lective modes (Bogoliubov-phonon mode). 
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FIG. 1: Polarization potential induced by fermion particle- 
hole excitation. Dashed line: boson (non-condensate), solid 
line: fermion. 



For the fermion Green's function in © , we use that for 
the unperturbed fermion single-particle state obtained 



from Eq. J3J : 

iG f (x u x 2 ) = ^V^i)^ 2 )*e-^ (tl - t2) 

3 

x [6{h - t 2 )6{e( - e F ) - 9(t 2 - t x )6{e F - ej)]. (9) 

To advance the calculation, we expand the boson field 
operators, (p and (p\ by boson creation/annihilation op- 
erators, P\ and pi for the quasi-particle states: 

X 
X 

where the quasi-particle wave functions, ux{z\) and 
vx{zi), for the eigenenergies h£l x can be determined from 
the Bogoliubov-type eigenequations: 



Cux(zi) + J dz 2 {n b g bb S(zi - z 2 ) + hH{z u z 2 ; Cl x )} {u x (z 2 ) - v x (z 2 )} = hn x u x {zi), (10) 
£-v x {zi) + J dz 2 {n b g bb S(z 1 - z 2 ) + 7ffl(z 1; z 2 ; Q x )} {v x (z 2 ) - u x (z 2 )} = -M^aO^i)) (11) 

I 



where C = K® + g bb n b + g b fUf. The hll denotes the 
energy-representation of the polarization potential in 
Eq. ©: 

P j h 

^t{zi)^p{zi)il)* p {z 2 )^ h {z 2 ) 



V'p (zi ) ^fe (zi ) (^2 ) Vfr 0*2 ) 



j j_ j 



(12) 



where the index p(h) represents a state with energy above 
(below) the fermi energy. Since we consider the station- 
ary condensates, the order parameter <!> is assumed to be 
real [28|. 

We should note that the number of eigenvalues fix ob- 
tained from Eqs. 1|1UI11|) may exceed the dimension of the 
bosonic quasi-particle space because these equations are 
not linear for the eigenvalue £l x (The Till depends on the 
eigenvalue). In this case, the eigenf unctions for different 
eigenvalues do not satisfy the orthogonality relations, but 
they are proved still to be a complete set. 

For normalizations of the quasi-particle wave func- 
tions, we take 

[dzi [\u x ( Zl )\ 2 - \vx( Zl )\ 2 } 



dz 1 dz 2 {u* x (zi) - v* x (zi)} 
dU(z l7 z 2 ; Q) 



dV. 



{u x (z 2 ) - v x {z 2 )} = 1. (13) 



III. 



PEIERLS INSTABILITY OF 
HOMOGENEOUS STATE 



A. Dispersion equation for Excitation energies 

Let's consider the case of u a = 0, and discuss the sta- 
bility of the ID homogeneous states in uniform mixtures, 
where the boson/fermion (ID) densities n b j and the or- 
der parameter $ take constant values: 



(14) 



where L is a quantization length for the box-potential 
regularization. In this case, the fermion single wave func- 
tions become plane waves, tpk = e lkz /^fh, with the wave 
number k = irn/L (n = integer) and the mean-field 



single-particle energy e{ = h k jlmj + g bb n b . The fermi 
wave- number and energy, k F and e F , becomes k F = nrif 
and e F = H 2 TT 2 n 2 j '(2m/) + g bb n b . 
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The boson quasi-particle wave functions u x and v x 
also become plane waves: u x ^ k (z) = u xk e lkz , v xk (z) — 
v\,ke- lkz ■ Substituting them into Eqs. (|1C)I11|> . we obtain 



(e\ - fJ, b )u\ tk + {gbbm + Kn. k (Q.\)} {u\,k - v\ tk } 



(15) 



Using the equality of the bosonic chemical potential and 
the interaction energy in the mean-field approximation, 
= gbbn-b + gb/rif, we obtain the dispersion equation for 



(hn x y = (elY + 2e b k {g bb n b + fin fe (fi A )}, (18) 



(el - Hb)v\ : k + {gbbm + hTl k (fi x )} {uA,fc - «A,fc} 



(16) 



where e k = H 2 k 2 /2m b + gbb^b + gbfUf, and the polariza- 
tion potential ftn^f^) for plane waves is given by 



hIl k (Q\) = gl f n b 



—kj 



dq 
2^ 



{ + irj 



fiO A - t{ + e f q+k + irj 



.(17) 



where e b k — h 2 k 2 /2m b . 

Introducing the complex energy £l x = tv x — we 
separate the polarization potential into the real and the 
imaginary parts 



hii k (u\n\) = nnf e (co x , lx ) + ihni m (io Xllx ). (19) 

After execution of the ^-integration in Eq. (|17|1 , they be- 
come ^3 



nn{ m (u; x , lx ) 



Ak 
2 



In 



In 



-A 



{{k + 2 k F )~ 2 ^} 2 + ( 2 -^) 2 { (fc + 2M + ^} 2 + (^) 

{(fc-2M + ^l 2 + (^) 2 
2mfL0 X 



{(fc-2fc F )-^> 2 + (^) 



n6 2k F 



arctan 



hk 

n / 2m f \ 9 

2 (-SF) W A7A 



n6 2k F 



2mfuj x 



hk 



arctan 



/ 2m t \ 9 
2 (-R^) W A7A 



(20) 



where Ai 



ifc = 9h ^. b , and the principal values should be 
taken for the arctangent function: — § < arctan a; < 
Using these results, Eq. (|18fl can be separated into the 
coupled equations 



h 2 ^ 2 --, 2 ) = (e b k ) 2 + 2e b k {g bb n b + hll^ 



WA7A 



4 



(21) 
(22) 



Now, in order to analyze the dispersion equations 
(|21I22|I qualitatively, we discriminate two cases. First, 
when ?ir4 m (w A ) j= and w x ,k > 0, Eqs. (J3TJ and |23 
have solutions with complex energies (7^ > 0), which 
correspond to the continuum p-h excitations in the con- 
tinuum. 

Second, in the region where any p-h pair excitations are 
prohibited (7^ = 0), the polarization potential becomes 
real: 



m? 1 (u>x) = -A k ]n 



(k + 2k F ) 2 - (2m f uj x /Tikf 



{k-2k F f - (2m f Lu x /hky 



(23) 



In this case, introducing the dimensionless variables uj x — 
2lu x / (v F k F ) and k = k/k F {v F — hk F /mf: the fermi 



velocity), the dispersion relation, eq. I|21|) . becomes 



~,2 



UJ 



m 2 ~ 4 
m b 



2k 



2 V B 



2- Sin 



k 2 (k + 2) 2 



k 2 (k - 2) 2 — uj 2 



(24) 

where vb = y/gbb^b/nib and £ is the dimensionless boson- 
fermion coupling constant defined by 



c 



9bf 



irgbbhvF 



(25) 



The solutions of Eq. (|24|) correspond to the bosonic ex- 
citation energies, which we are interested in, and vb is 
just the Bogoliubov-phonon sound velocity. It should be 
noted that the scaled equation depends on three 

dimensionless parameters, irif/mb, vb/v f , and £. 



B. Excitation Spectra and Peierls Instabilities 

Let's study the bosonic excitation spectra on the homo- 
geneous ground state (|14|) with solving the dispersion re- 
lations H21I22(I . To understand general features, we show 
numerical results for two typical cases, (a) vb < v F and 
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(a) v B < v F 




(b) v B > v F 




0.0 0.5 



FIG. 2: Excitation spectrum (in units of fefiiF) of bose- 
fermi mixtures. Solid lines: collective modes of p-h pairs and 
bosonic modes, dashed line: free boson, shaded area: contin- 
uum p-h pair excitation spectra. 



(b) v B > v F , in Figure. 2: (a) (v B /v F ) 2 = 0.67, C = 0.99 
and (b) (vs/vp) 2 — 1.50, C = 0.80. The same values 
have been taken for boson and fermion masses, m b = m/, 
for all cases. The shaded areas in Figure 2 correspond 
to the continuum spectra corresponding to the fermionic 
p-h excitation states, and the solid lines are to the iso- 
lated modes, corresponding to the collective excitations. 
It should be noted that two collective modes (low- and 
high-lying) exist with the same wave number k (except 
narrow regions around 2kp)', they can be interpreted 
to be coherent superposition of the Bogoliubov-phonon 
mode from bosonic atoms and the fermion p-h collective 
mode. For comparison, the Bogoliubov-phonon spectra 
in boson- fermion nonintcracting systems (g b f = 0) have 
been plotted in Figure 2, and we can find that it runs 
between two collective spectra in the interacting ones. 

Let's discuss the specific features of the collective 
modes, especially in small k and k ~ 2kp regions. 

In small wave-number region, k < kp\l — (vb /vf) 2 \, 
the Bogoliubov-phonon modes (dashed line) in the 
boson-fermion noninteracting limit are located below the 
fermion p-h excitation modes in energy and absorbed into 
the fermion p-h continuum when vb < vp (case (a)). In 



case that vb > vp (case (b)), the corresponding modes 
are above the continuum. It suggests that the low-lying 
collective mode in (a) is mainly a Bogoliubov-phonon 
mode from bosonic atoms but the corresponding mode 
in (b) is mainly the fermion p-h one, and the high-lying 
modes have opposite characters. In Figure 2, the bosonic 
excitation modes are found to be influenced by the inter- 
action with the fermion p-h excitation modes, and pushed 
downward/upward in case (a) /(b). 

Around k ~ 2kp, the low-lying collective modes be- 
come very soft and, between two critical wave-numbers 
fc_ < k < k + , they are found to be zero-modes in both 
cases. To calculate the energy of these modes analyti- 
cally, we expand Eq. (|22|l to the leading order of \q\/kp 
{q = k- 2k F ): 



h 2 oj\~f\ = A 2 k F e\ [ir6{u\ - v F \q\) 



— arctan 



vpq 2 



(26) 



When lu\ < vp\q\, the above equation reduces to wa7a = 
0. In the case of u \ = "fx = 0, two critical wave-numbers 
are obtained from Eq. 121|) : 



fc± 
kp 



2±4exp 



C 



1 



rrifVp 
m b v B 



(27) 



which gives (fc_ = 1.97fci?, k + = 2.02kp) for (a) in Fig- 
ure 2, and (k- = 1.92fc F , k + = 2.05k F ) for (b) in the 
same figure. 

Analyzing the eigenvalues of Eq. (|21J) more detail, we 
find that they become purely imaginary (lu\ = 0, 7a ^ 
0) in the region between the two critical wave-numbers, 
< k < k + ; it is just the "Peierls instability" in Q1D 
system, which suggests that the homogeneous state is 
unstable against the fluctuations with the wave-number 
around 2k p. 

Finally, we should comment on a different kind of in- 
stability that may occur in the small wave-number fluc- 
tuation jig . When we take the k — > limit with keeping 
uj\/k — constant, the real part of the polarization poten- 
tial becomes 



an£> A ) 



(g bf ) 2 n b {vFk) 1 



Trhvp 



{v F kf 



(28) 



and the imaginary part n{ m vanishes except at vb = 
vf. Substituting Eq. 128fl into Eq. I|21|l . we obtain two 
branches of excitations: 



Hlu± = fik 



(7 



„2 \2 



) 2 + 4(v 



F V B 



1/2 



(29) 

In the strong boson-fermion interacting case (C > 1), the 
energy of the low-lying mode becomes purely imaginary, 
and it suggests an instability of the system. We find that 
the stability condition for the k ~ mode: 



-nt > 



(9bfl 
9bb 



(30) 
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where we have used the relation rif = hp /it. This k ~ 
instability is considered to cause the phase-separated 
states more stable. In the present paper, we do not dis- 
cuss this instability and concentrate only on the density- 
wave states caused by the Peierls instability, which is 
discussed in the next section. It should be a very in- 
teresting problem to study the cooperation/competition 
between these states. 



IV. DENSITY WAVES IN BOSON-FERMION 
GROUND STATE 

The results obtained in the previous section shows that 
the homogeneous state is not the true ground state of the 
system and a lower energy can be obtained for the state 
with a spatially-periodic condensate characterized by the 
2fci?-periodic order parameter: $(2) = b u + b p cos (2k fz). 
The constant bosonic density is given by m — b 2 + bp/2. 
In the case of weakly-periodic variation b p <§C b u , up to 
second order in b p , the mean-field Hamiltonian for the 



bose-fermi system becomes 



H 
T 



9bf 



where Ck and c k are the annihilation and creation opera- 
tors for fermions with the wave-number k and the energy 
e{ = h 2 k 2 /2m f + g bf (bl + b 2 p /2). The e b p in Eq. (23) 
is the energy generated from the periodic condensate: 
e b p = e b 2kp /2 + gbbb 2 , . The last term in Eq. (J31J is for the 
processes that the fermion with a wave-number k feeling 
periodic- and uniform-bosonic condensates scatters into 
the sates with k ± 2k p. 

In order to obtain a single-particle energy, we calculate 
the fermion Green's function: 

iFk, k >(t-t') = {T[c k (t)cUt')}) = J ^V-C'-Oi^), 

(32) 

where c k (t) = e lHt / h c k e~' lHt / h . For a fixed wave-number 
k, the Green's function F ktk i (u>) has non- vanishing off- 
diagonal matrix elements at F kjk ±2 kF in the Hamiltonian 
(|31|l . In the limit k — + +kp, fermion states with wave- 
numbers k and k — 2kp are almost degenerate in energy, 
but that with k + 2kp is separated. Consequently, in 
case of kp > k > 0, the Dyson equation for the Green's 
function can be approximated by 



!) » i pE( £ i-2fcF^ +4+2fc F £ fc). ( 31 ) 



hu) 



J 



-A 



-A Tiu — e 



k—2kr. 




(33) 

where F ktk is a diagonal part of the Green's function, and 
A = gbfb u b p . They can be solved easily: 



Fk,k — 



LU-E k /h—ii] uj — E k /h + ir] 



,(34) 



F k , k -2k t 

where 



-u k v k 



UkV k 



uj — E, /h — ir/ to — E k /h + if] 



U k = 



Vu 




(35) 



(36) 



(37) 



with £fc = hvp(k — kp). The single-particle energy with 
wave- number k is obtained from the pole of i*k,fc(w): 



E. 



(38) 



from which we can understand that the 2\ A| is the energy 
gap at the fermi surface. It is clear that the periodic 
condensate has produced an energy gap 2|A| near the 
fermi surface. 

To evaluate the value of A, we calculate the total en- 
ergy of the system in the present approximation. The 
fermionic contributions to the energy is obtained by 
summing-up the E k in Eq. I)38fl up to the fermi wave- 
number kp: 



E, = j dkE k 6(kp - \k\) = ( y + 9b f n b ) n f 



+ n f 



e F 




+ v / 4 + A5 



(39) 



where eF — hvpkp/2. For simplicity, we consider the 
weak coupling limit, A -C e/r, then the total energy den- 
sity of the system becomes 



E tot (\A\) ~ E tot (0) 



4e? 



C 



rrifVF 
m b v B 



1 f 4e F 

2 r +2in w 



(40) 



Minimizing Eq. (|40|l with respect to |A|, we obtain the 
stationary value of the gap parameter 



^=4exp 
eF 




(41) 



Using the above result, the energy difference between the 
ground and the uniform (A = 0) states is 



E tot (\A\) - E tot (0) = - 



A 2 



(42) 



which shows that the state with the "periodic conden- 
sate" (bp 7^ 0) has lower energy than the uniform one. 
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FIG. 3: Gap energy as function of rrijVp jm^vs (in units of 
e F ) for C = 0.99. 



Next, we discuss the fermionic density, which is ob- 
tained from the Green's function: 



+Fk-2k F ,k£ 



-i2kf z 



Fk,k + Fk,k~2k F e 



+ Fk-2k F ,k-2k F 



i2k F z 



(43) 



Using Eqs. I|34I35I36I37|I . it becomes 

H f( Z ) = Y} U k + V k - 2U k V k COS (2k F z)} 
k 

A 



= n f 



1 



cos (2kpz) 



(44) 



which shows that, for finite A, the fcrmion state is in 
the density-wave one where the fermionic density has a 
periodicity 2kp. Corresponding to Eq. I|44|) . the direct 
calculation shows that the bosonic density rib(z) has the 
similar density-wave structure: 



n b (z) 



in, 



\9bf\n b 



cos(2fcpz) 



(45) 



It should be noted that, different from n/(z), the peri- 
odic part in rib(z) is proportional to the absolute value of 
the gap parameter |A|. That means that the boson and 
fermion density-waves are out of phase for the boson- 
fcrmion repulsive system, but in phase for the attractive 
one. 

As seen in Eq. 1|41[1 . the gap energy A essentially de- 
pends on two combinations of parameters, rrifVF /rribVB 
and the dimensionless coupling constant £. In Fig- 
ure. 01 we just plotted the gap energy as a function of 
rrifVF /mbVB for £ = 0.99, from which we can read off 
that the boson-fermion density-wave state can be ob- 
served in increasing the ID density of condensate and/or 
using heavier bosons than fermions. 



V. PEIERLS INSTABILITY OF SYSTEM IN 
HARMONIC OSCILLATOR POTENTIAL 

In this section, we briefly discuss the non-uniform 
boson-fermion mixtures confined in a HO potential with 
the finite axial dimension (u a ^ 0), and the condition for 
its Peierls instability. 

Using the eigenenergy = hui a {n + h) (n = 

0, 1, 2, • • •) and the wave function (f>^° for the ID HO po- 
tential, the quasiparticle amplitudes u\ and v\ are ex- 
panded by 

UA (z) = 5>^°(z), v x ( Z ) = 5>^°(2)- (46) 



Substituting them into Eqs. (|10lll|l . we can obtain the 
eigenequations for the excitation energies ftwy. 

Y, W h n° - riujx)5 m , n u x n + nn mn {^ x ){u x n - v x n )] = 0,(47) 

n 

Wn + hu>x)5 m , n v* + hIL mn (ux)(vtl - u*)] = 0(48) 



where, to concentrate on the role of the polarization po- 
tential effects in finite systems, we have neglected the 
Hartree potential. The polarization potential hIL mn in 
Eqs. (|47I48[) . is given by 



hU mn (ujx) = {gbf) 2 N b 

mh\p0)(0p\hn) (mp\hO) (Oh\pn) 



><E 

ph 



hujx 



The matrix elements in Eq. I|49|l are defined by 



ho 



(mh\pO) = / dz^CV£V°> 



(49) 



(50) 



where the single-particle wave functions in Eq. (|50|l have 
been replaced by the HO eigenfunctions (fif . It should be 
noted that the matrix element Eq. 1)50(1 has a transition 
probability only in the spatial region where condensate 
exists. It is important at this point to contrast finite 
system to homogeneous system in which condensate is 
present all over the range. 

We restrict the present discussion into the parameter 
region: %uj r > pf > fib > TioJa- Under these conditions, 
we can take several approximations. Because of fib > 
huj a , the system is in the Thomas- Fermi (TF) regime, so 
that we can take the TF approximation for the order- 
parameter: 



$ = ,/Nb<poocJl-- s -0(zTF-\z\). (51) 



TF 



The ztf in Eq. (|51|l is a cut-off length in the TF approx- 
imation: 



z T f 



(3iV [3\ 



(52) 
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where = mg^aho / 'h 2 with aho = \Z^/ mu « (the HO 
length in the axial direction). 

As we have seen in chapter III, Peierls instability 
comes from the coupling between the boson mode with 
wave-number near 2kp and the fermion p-h excitations 
near kp . Although the wave- number is not strictly con- 
served in a HO potential, the instability mode would have 
'wavenumber': k n = \/2n/aho ~ 2kp. Since the higher 
nodal modes n ~ Nf have a broader spatial extension 
(L n = \pZnaho) than a condensate one, the asymptotic 
expansions can be used for the </>^° in the middle of the 
trap \z\ < ztf < L n |3Cj : 



t%°(z) (x cos (k n z — 



(53) 



For a bosonic eigenstate with higher nodal mode, the 
negative energy amplitude = of Eq. (|T7|) can be 
neglected. 

Using these approximations, the matrix elements 
like (Op\hn) become linear combinations of the form 
Ji(ztfK)/ztfK, where J\(x) is the lst-order Bessel 
function and the if is a wave-number difference between 
the initial and the finial states: K — k n + kh — k p . In 
the case of |x| 3> 1, using the asymptotic expansion 
of the Bessel function, J\(x)/x ~ | | 3/2 ^ terms f 
\x\ < 0(1) can contribute to the matrix elements. Under 
the restriction of wave- numbers kh < kp and k p > kp 
with k n ~ 2kp (kp = w2Nf/ah ), the main contribu- 
tions of matrix element are from the scattering processes: 



kn. 



kh and k n + kh — > k p . To make an estimate, 



we replace a product of Bessel functions as 



Ji(xi) Ji(x 2 ) 



Xl 



X 2 



4 V 



\ Xl \)0(l - \ Xl - x 2 \), (54) 



and take the continuum-limit for the wave-number sums: 

Finally, the eigenenergy for a collective state with the 
wave number k n = 2kp can be obtained 



Tiuj\ = 4ei 



gb/^b(0) 
2-kTivf 



ln |4/cf-Ztf| 



(55) 



Contrary to the uniform system, the static polarization 
potential, the second term of right hand side of Eq. (1551) , 
is not divergent at 2k p. This is because the scatter- 
ing of trapped atoms without wave-number conservation 
smears a singularity due to a sharp fermi sea. The Peierls 
instability occurs where the polarization energy is over- 
come by the kinetic energy Aep: 



1 < 



<>b(0) 2 



ln |4fcp ztf | 



(56) 



where vb(0) — y/gbbnb(0)/m. This is the Peierls insta- 
bility condition where boson-fermion density wave may 
occur. 

We roughly estimate the parameters according to ex- 
perimental conditions. When mb/mp ~ 1, we can ex- 
press the velocity ratio and the dimensionless coupling 
constant by 



AN, 



3N b 



c = 



u r abb 
u a aho 



n 2 Nf uj a abbdho 



(57) 
(58) 



As a possible candidate for the Peierls instability, we 
take the rubidium isotope system: 87 Rb— 84 Rb mixtures. 
Taking the scattering lengths in 3lj, abb — 5.3 nm and 
dbf = 29.1 nm and the typical trapping frequencies, 
uj a = 2tt x 10 Hz and uj r — 2tt x 15 kHz, Eq. (|56*|) gives 



N 



f 



W and N h 



10 for the realization of the 



Peierls instability. 



VI. SUMMARY 

In the present paper, we studied the occurrence 
of Peierls instability in Q1D bose-fermi mixtures at 
zero temperature. We analyzed the bosonic collective- 
excitation spectra in random phase approximation. It 
shows that the mixtures of uniform BEC and a fermi gas 
are unstable against a spontaneous formation of a collec- 
tive mode of wave-number 2k p; this type of instability 
is known as Peierls instability. This result suggests that 
the ground state of bosons is a periodic condensate with 
the period ir/kp. In the variational method, the boson- 
fermion density wave state have been shown to have a 
lower energy than the uniform state. We also expanded 
our analysis for systems in an axial harmonic oscillator 
potential and derived Peierls instability condition. 

It is well known that density waves in Q1D conduc- 
tors lead to generations of a Lee-Rice- Anderson |32| and 
a phase soliton _24] modes. Thus, studies of the dynam- 
ical properties of boson-fermion density waves should be 
interesting problems in future. 

The authors are supported by Grant-in- Aids for Scien- 
tific Research by JSPAS. One of the authors (Miyakawa) 
was also supported by JSPS fellowship of Japan. 
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